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ABSTRACT 



C . I 



A discrete-event simulation of a stochastic process, a 
machine-repairman model, has been programmed on the IBM 
Personal Computer. The model consists of three helicop- 
ters, of which two are in service and one is in cold 
standby, with an option of one or two repairmen. 

The program output is a graphics display containing a 
system state-versus-time graph, a table of statistics, and 
animated figures that illustrate the current state of the 
system. The program user can directly observe the dynamics 
of the model as the fixed-increment, simulation clock 
advances. The user has the option of changing the follow- 
ing model parameters: helicopter failure rate, repairman 

service rate, and the number of repairmen to employ. 
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I. 



INTRODUCTION 



During recent years, the versatile microcomputer has 
found a niche in many of the professional discliplines 
including education. For difficult and complex subjects, 
educators are exploring the possibilities of the microcom- 
puter as a teaching aid in addition to its more traditional 
role as a computing machine. Some universities now use 
microcomputer graphics to illustrate principles that cannot 
be fully explained in a lecture. [Ref. 1] 

A widely used application of computers is to imitate or 
simulate stochastic processes. This technique usually 
involves only the numerical evaluation over time of the 
corresponding model with the intent of gaining some under- 
standing of the process. But using graphics as part of the 
output in a simulation provides the user another 
perspective. For example, observing a working model--such 
as a birth-death process — can help in grasping the funda- 
mentals of the model and provide some insights into the 
process not easily seen on the classroom chalkboard. 

This paper describes the development of a discrete- 
event, graphic simulation of a specific example of the 
machine-repairman model. The simulation gives the observer 
the impression of real-time passage. A graph on the visual 
display shows the history of the system's transitions from 
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state to state, and animated figures illustrate the 
dynamics of the model. The simulation also exhibits a 
table of the theoretical, and the current, estimated values 
of quantities commonly used in evaluating queueing systems. 

Section II of this paper is a discussion of the 
machine-repairman model, the assumptions upon which the 
model is based, and the methods used - for obtaining the 
theoretical and estimated values of the quantities of 
interest. Section III covers the simulation's implementa- 
tion on the IBM Personal Computer (IBM PC). Appendix A is 
a user's guide which contains summary information about the 
machine-repairman model and instructions for operating the 
program. The program flowchart and the program listings 
are found in Appendixes B, C, D, and E. 
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II. THE MACHINE-REPAIRMAN MODELS 



A. DISCUSSION 

The simulation described in this paper actually con- 
tains two machine-repairman models. In both models .there 
are three machines, represented • by helicopters. The models 
are distinguishable by the number of repairmen; there is 
one repairman in the first and two repairmen in the second. 

The nominal mission is to keep two helicopters flying 
at all times. The third helicopter is placed in cold 
standby as a backup to the other two. If available, a 
repairman is assigned immediately to a helicopter upon its 
failure. If a repairman is not free, then the helicopter 
joins a queue to await repair. 

Several assumptions are made to simplify the models. In 
the two-repairmen model, repairmen are assumed to be 
equally competent and to take the same mean time to repair 
a helicopter. Both models assume the repair time to be 
continuous; that is, repair parts are always available, so 
there is no delay in repairing a machine other than the 
time required by the repairman to perform the work. The 
repairmen work independently and do not assist one another. 
Similarly, the models assume the helicopters operate inde- 
pendently and have the same service life with the same mean 
time to failure. By assumption, a helicopter is not 
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subject to failure while in cold standby since it is not 
operating . 

A machine-repairman model typically assumes its 
machines operate continuously until failure. In the case 
of aircraft, this assumption does not reflect reality since 
in actual practice aircraft do not fly continuously,’ but 
have instead periods of inactivity between missions of 
limited duration. As already noted, one can reasonably 
assume that an inactive helicopter is not subject to fail- 
ure and that the failure process occurs only when the 
helicopter is active. 

The actual failure rate of a flying helicopter is the 
number of failures per flight hour. The product of the 
actual failure rate times the ratio of the number of flight 
hours performed to the number of calendar hours elapsed 
until failure defines the helicopter's effective failure 
rate. Thus the effective failure rate is the number of 
failures per calendar hour. To simplify the simulation, 
the models assume that the helicopters operate continuously 
until failure with the effective failure rate. 

The selection of helicopters to represent the machines 
in these models allows a simple graphic design in depicting 
their status. An operating machine is seen in the display 
as a green, flying helicopter; a helicopter that is in cold 
standby is also colored green but is not flying. A machine 
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that is down and awaiting repair is colored red. When the 
down machine is in repair, it is seen with a repairman. 

Figure 1 is a frequently used, pictorial representation 
(see [Ref. 2], for instance) of the machine-repairman model 
with one repairman. The dashed lines enclose the model's 
repair-queue system which contains a single-server queue. 
The machines outside the dashed box are the potential 
customers for the repair queue. Since there are only three 
machines in the model, the longest possible length of the 
queue will be two. 




Figure 1: Machine-Repairman Model with One Repairman. 

This model has a second queue, the cold-standby queue, 
since a helicopter may have to wait before being assigned 
to a mission. The standby queue can contain only one heli- 
copter, and this occurs when there are no machines in the 
repair-queue system. 
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The cyclic nature of the model is clearly evident in 
the figure. A helicopter's life cycle consists of an oper- 
ating period followed by, after a possible delay in queue, 
a repair period which is in turn followed by, after another 
possible delay in cold standby, another operating period. 




Figure 2: Machine-Repairman Model with Two Repairmen 



Figure 2 depicts a pictorial representation of the 
machine-repairman model with two repairmen [Ref. 3]. In 
this system, helicopters waiting in the repair queue go to 
the first available repairman for service. The maximum 
possible length of the repair queue is only one. As in the 
case of the one-repairman model, the standby queue can 
contain at most one helicopter. 

The machine- repa i rman model is a special case of a 
birth-death process in which the event of a machine failing 
at time t constitutes a birth. Let T be the operating life 
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of a helicopter until failure and assume that T is exponen- 
tially distributed with a failure rate a. Then the mean 
time to failure (MTTF) of a helicopter is E(T) = 1/a. 
Similarly, the completion of repair on a machine is a 
death, so define S to be the service time of a repairman. 



The service 


time is 


also 


assumed to 


have 


an exponential 


distribution 


with a 


service rate p. 


The 


mean time for 


repair (MTFR) 


is then 


E(S) 


= 1/u. 







For these machine-repairman models, the state of the 
system, N(t) = n, is defined as the number of helicopters 
down at time t. Since there are only three helicopters in 
the models, then n = 0,1, 2, 3. 



B. MODEL ESTIMATORS 

In the study of queueing models, there are four quanti- 
ties that are commonly used to evaluate the queueing sys- 
tem's performance. As part of the graphical output, the 
simulation will provide the user a table of theoretical and 
estimated values of these four quantities: 



L 

L 

q 

W 

w 

q 



the average 
system , 


number 


of 


machines 


in 


the 


repai r 


the 


average 


number 


of 


machines 


in 


the 


repair 



queue , 

the average amount of time a machine is down, 

the average amount of time a machine waits for 
repair . 
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As shown in Figures 1 and 2, the variable q is the time 



that a down machine 
W = E(q). It follows 

q 

helicopter is in the 
average waiting time 
time : W = W + E (S ) . 

q 

Little's formulas 
quantities of interest 



waits in the repair queue. Then 
that the average total time that a 
repair system is the sum of the 
in queue and the expected service 

show the relationship of the four 



L = AW 

L rr = 

q q 

where X is defined as the average arrival rate of machines 
entering the repair-queue system. 

There are three additional quantities that will be 
computed and presented by the simulation. They are defined 
as 



= P {down machine must wait for repair}, 

P g = P{up machine must go to cold standby}, 

Av = availability = the proportion of time that at 
least one helicopter is not down. 

To calculate all of eight of these values, it is first 
necessary to compute the system's steady-state probabili- 
ties, p n =P{N=n}, n = 0,1, 2, 3. In other words, p n is 
the long-run probability that the system is in state n. The 
details of calculating p^ , a function of the machine 
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failure rate and the repairman service rate, will not be 
discussed in this paper. The interested reader should 
refer to the bibliography for a list of textbooks on queue- 
ing theory and birth-death processes. In addition to the 
limiting probabilities, {p n , n = 0,1,2, 3}, the program 
computes the fraction of the time the system is in state n, 
denoted p^ , as equal to the total time in state n divided 
by the total elapsed time. The quantity p n is an estimate 
of p n . 

The value of L, the theoretical average number of 
machines in the repair system, is found by averaging the 
number of down machines over all states: 



L = ) n P n = (UP-. + (2)p ? + (3)p . 

n^O n 1 Z J 



The estimate of L, denoted L, is computed with the same 
expression with the exception that the estimate p^ is sub- 
stituted for p . 

^n 

The average number of machines in the repair queue L q 
and its estimate L q depends in part upon the number of 
repairmen in the model. In the one-repairman case, 
machines occupy the repair queue only when the system is in 
state 2 or 3 , thus 



L q = P 2 + 2p 3 , and 
L q = f 2 * 2^3 • 
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In the two-repairmen case, a queue forms only when the 
system is in state 3, thus 



L 



q 



L 



q 




and 



The average rate at which helicopters enter the repair 
queue, X, is dependent upon the number of active machines 
subject to failure in state n and the proportion of time 
the system is in state n. There are two active helicopters 
in states 0 and 1, one in state 2, and none in state 3. The 
combined failure rate of the active helicopters in state n 
is given by 



2 a n = 0 , 1 



a n = a n = 2 



0 n = 3 



Then X is found by 



X = ] oi p„ 
n = 0 n n 



= 2ap 0 + 2ap 1 + ap 2 



For the estimate X, the program divides the number of fail- 
ures that have occurred in time t by t , that is 



X = total number of failures/total elapsed time. 
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The average waiting time quantities, W and are then 
calculated by using Little's formulas 



W = L/X 





the length of time each helicopter spends in the repair 
queue and uses the expression 



where q^ is the waiting time in queue of the ith helicopter 
and m is the total number of helicopters that have entered 
the queue by time t. Then W is easily found with 



The estimate of the probability that a down machine 
must wait for repair is the proportion of the time that the 
system spends in a state in which a machine is down with no 
unemployed repairman. The repairman is not available in 
the one-repairman model when the system is in state 2 or 3 : 



In the two-repairmen case, the repairmen are busy when the 
system is in state 3: 



,m 



W 



q 



(1/m) ) q . m = 0 , 1 , 2 

i = l 1 



J • • • 



W = W + E(S ) . 

q 



p 



d 
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For P , the only time a machine must wait in cold standby 
is when the system is in state 0; therefore, P g = p^. 

The availibility of at least one helicopter in the 
system is the sum over all states n of the percentage of 
the flyable helicopters in each state multiplied by the 
proportion of time the system is in that state; 

Av = Pq( 1) + P]_(2/3) + P2<l/3) + p^tO). 

The estimate of the availabiliy is calculated by substitut- 
ing p n , n = 0,1, 2, 3, into the above expression. 
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III. THE SIMULATION 



A. THE COMPUTER 

The microcomputer used for this simulation was the IBM 
Personal Computer. Since the program MACHREPR uses color 
graphics extensively, the computer requires a color display 
monitor and either IBM's Color/Graphics Adapter or a suit- 
able graphics card from another manufacturer. The program 
will not run without a graphics adapter. Other minimum, 
hardware requirements for this simulation are 128 kilobytes 
of read-write memory (RAM) and at least one 5i-inch, floppy 
disk-drive . 

At the beginning of this project, the author chose to 
use IBM's disk operating system PC-DOS with the intent of 
enabling the simulation to run on other microcomputers 
using the more generic version of PC-DOS: MS-DOS by 

Microsoft. It soon became apparent, however, that MACHREPR 
would need to directly access certain parts of the 
computer’s memory in order to improve the program's 
performance. Since the program makes direct calls to the 
PC-DOS, there can be no assurance that the simulation will 
run on other "IBM-compatible" machines. 
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B. THE PROGRAMMING LANGUAGES 



MACHREPR is written in Advanced BASIC which is the pro- 
gramming language supplied with PC-DOS. Advanced BASIC was 
selected because of its convenient graphics commands and 
its ready availability. Although Advanced BASIC is an 
interpreter, program speed was not an essential criterion 
that warranted using a compiled language. Advanced BASIC 
was insufficient, however, in two aspects, and so MACHREPR 
uses two subroutines written in 8088 assembly language to 
improve program performance. These two subroutines, named 
SCRNSHFT and RNGEN, will be described later in this 
section . 

C. PROGRAM DESCRIPTION 

Appendix B contains the flowcharts for the BASIC pro- 
gram MACHREPR. The flowcharts illustrate the simulation's 
general construction and program flow; however, they are 
not a verbatim duplication of the program listing. The 
program listing for MACHREPR and the two source-code list- 
ings for RNGEN and SCRNSHFT are in Appendixes C, D, and E 
respectively . 

After initializing several variables and arrays, the 
program presents the title screen. The program menu 

appears next with the following options: (1) see the pro- 

gram instructions, (2) change the model's parameters, (3) 
use the default parameters, (4) set the random-number 
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generator seed, or (5) end the program. The instructions 
that are available on the screen are a summarized version 
of those found in Appendix A, the user's guide. The model 
parameters that can be changed by the user are the mean 
time to failure, the mean time for repair, and the number 
of repairmen. 




Figure 3: An Example of the Graphics Display. 



Figure 3 shows an example of the graphics display 
screen as it appears during the simulation. The upper half 
of the screen contains the system state-versus-time graph. 
The system state values are on the ordinate of the graph, 
and fifty divisions representing time are on the abscissa. 
The three graphic figures of helicopters in the lower left 
corner symbolize the status of the system. The helicopter 
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seen above the others is in active service; the two 
machines with broken rotors are down. In this example, the 
helicopter graphics indicate that the system is in state 2, 
where helicopter B is flying, helicopter C is under repair, 
and helicopter A is awaiting repair. The graph shows that 
the current time is 69 and that the system jumped from 
state 1 to 2 at t = 68. 

The table in the lower right corner contains the esti- 
mated and the limiting values of the quantities discussed 
in Section II. The row labels are: 

Ar - the arrival rate of machines into the repair 

queue , 

Wq - the average wait time of a machine in the repair 
queue , 

Lq - the average number of machines in the repair 

. queue, 

W - the average amount of time a machine is down, 

L - the average number of machines in the repair 

system, 

Avl - the proportion of time that at least one machine 
is available. 

Pd - the probability that a down machine must wait for 
repair , 

Ps - the probability that an up machine must go to 
cold standby. 

MACHREPR consists of four major parts: (1) main pro- 

gram, (2) clock module, (3) failure module, and (4) repair 
module. The main program handles the housekeeping details 
such as initializing variables, loading the machine 
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language subroutines, generating the graphics display, and 
performing the tasks selected from the program menu. Within 
the main program are routines for computing the limiting 
values of the tabulated quantities and for printing them on 
the display screen. The main program also calculates the 
initial failure times for the first entry into the clock 
module . 

The clock module determines the next-event time, 
whether the event is a failure or the completion of repair, 
and which helicopter it affects. The simulation clock is a 
fixed-increment, time-advance type. In contrast to a next- 
event, time-advance clock which skips over the intervening 
time between events, the clock module in MACHREPR advances 
the time in equal increments. With this type of clock, the 
user experiences the illusion of real-time passage. 

During the execution of the simulation, the clock 
module monitors the keyboard for certain commands issued by 
the user. Pressing the "P" key pauses the program to allow 
the user to study the graphic display. Pressing the "C" 
key causes the program to continue. The "S" key will stop 
the simulation and return the program to the program menu. 
Should the clock reach the preset limit of 9950 time units, 
the program will automatically stop the simulation and 
return to the program menu. 

Every time the clock is advanced one increment, the 
clock module updates the system state-versus-time graph 
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according to the current state and time values. As the 
time advances, the plotted lines will eventually reach the 
edge of the graph. When this condition is met, the clock 
module shifts the entire graph left one time unit and draws 
the new time increment. This procedure continues for all 
successive time increments. When the simulation clock has 
reached the time for the next event, the program jumps to 
either the failure module or the repair module as 
appropriate . 

The failure module changes the helicopter display 
according to the current state and model parameters . The 
program keeps track of the machines individually, so that 
the graphic image of the particular, failing helicopter 
transforms from a green, flying figure to a red, grounded 
one. If a repairman is free, the graphic figure of a man 
is drawn next to the aircraft. If no repairman is avail- 
able, then the human figure is omitted, signifying that the 
helicopter is placed in the repair queue. If a repairman 
was assigned, the failure module calls the random-number 
generator and computes the service time. The sum of the 
service time and the current time is the time that the 
repair will be completed and this is compared to the next- 
event times of the other machines after the program returns 
to the clock module. 

In a similiar manner, the repair module changes the 
graphic figure of the helicopter that has just completed 
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repairs. If there are less than two helicopters flying 
when the repair event occurs, then the repaired machine 
transforms from red colored and disabled to green colored 
and flying. If two helicopters are already flying, then 
the repaired machine goes into the standby queue. The now 
free repairman moves to the helicopter at the head of the 
queue, or if the queue is empty, he disappears from the 
screen . 

Before the program returns to the clock module, the 
repair module computes the interarrival time to failure of 
the helicopter just repaired. Failure times are not com- 
puted for helicopters placed in the standby queue. If the 
repairman is assigned to another helicopter, the repair 
module also computes the service time. 

After returning to the clock module, the simulation 
updates the estimates of the tabulated quantities and 
prints them on the display screen. Both the failure and 
repair modules change the state-indicator variable to its 
new value so that the clock module will draw a horizontal 
line one time unit in length at the proper place on the 
graph . 

SCRNSHFT is the subroutine that shifts the system 
state-versus-time graph one time unit to the left per clock 
cycle. The subroutine uses certain program transfers 
called interrupts to the PC-DOS Basic Input/Output System 
(BIOS) which contains routines that control the video 
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display. SCRNSHFT is not appreciably faster than a BASIC 
routine designed to perform the same task; however, the 
subroutine eliminates an annoying flicker in the display 
that was characteristic of the BASIC version. 

D. PROGRAMMING CONSIDERATIONS 

The program computes the interarrival (failure) and 
service times, assumed to be continuous, exponentially dis- 
tributed, random variables, in the single-precision format. 
The current time, an integer value, is added to the inter- 
arrival time and to the service time to yield the next- 
failure time and the next-repair time respectively. In the 
unlikely event that the two single-precision, random 
variables match, the simulation will proceed as if the 
failure occurred earlier. 

Whenever the model jumps from one state to another, the 
program computes a next-failure time and a next-repair time 
for the helicopter just affected by the transition. Thus 
MACHREPR always keeps in memory six next-event times for 
all machines regardless of their individual status. To 
determine the next event, the clock module compares the 
smallest of the three failure times with the smallest of 
the three repair times. With this method, there must be a 
means of preventing the program from simulating events con- 
sidered impossible by the assumptions of the model. 
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For example, one model assumption is that a helicopter 
in cold standby is not subject to failure. It follows that 
the same machine cannot be repaired since it is not in the 
repair queue. Similarly, a helicopter that is in repair 
cannot fail and a helicopter that is flying cannot be 

repaired. To handle these impossible events, the program 
will set the failure time or the repair time, as appropri- 
ate, to equal 10^ which represents infinity. 

Because of the discrete nature of the graphic display, 
all time values must be converted to integers prior to 
plotting the state-versus- time graph. It is for this rea- 
son that the random number representing the next-event 
time, be it a failure or a repair, is rounded up to the 

next higher integer. The result is that the interarrival 

and service times are not exponentially distributed but. are 
instead geometrically distributed. 

The discreteness of the time variable introduces a bias 
in the estimators since they are based on the assumption of 
continuous distributions. The effects of the bias, how- 
ever, are minimal. The values of the estimates are still 
seen to converge to the limiting values as time advances. 
Despite the generation of geometric random variables, the 

simulation still provides a sufficient visual demonstration 
of the machine-repairman model. 
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E. RANDOM NUMBER GENERATION 

The author found BASIC's random-number generator, a 
function called RND , unsuitable for simulation. Even so, a 
sample of 10,000 uniform variates produced by RND passed 
with .95 confidence three standard, nonparametr ic tests for 
uniformity and independence: the serial test, the frequency 
test, and the runs up-and-down test. On the other hand, a 
two-dimensional scatter plot of 5000 pairs of uniform vari- 
ates, shown in Figure 4, clearly illustrate the lattice 
structure of the RND generator. 




Figure 4: Two-Dimensional Plot of 5000 Pairs of Uniform 
Variates from the BASIC RND Function. 

MACHREPR uses a random-number generator called RNGEN . 
Written in 8088 assembly language, RNGEN is modification of 
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a generator programmed by Associate Professor B.O. Shubert 
of the Naval Postgraduate School for use on the IBM PC. 
RNGEN uses the algorithm described by Fishman [Ref. 4] 
which is based on the expression 

X. = 16807*X. mod (2 31 - 1) 
l+l l 

where X. is the old seed and X. , is the new seed, 
l l+l 

Dividing the new seed by the modulus produces the uniform 
(0,1), random number. 

This algorithm is the one used by LLRANDOM, a random- 
number generator developed at the Naval Postgraduate School 
and described by Lewis, Goodman, and Miller [Ref. 5] as 
able to produce good quality, uniform variates. 

The simulation computes exponential random numbers by 
the inverse probability integral transformation of uniform 
variates produced by RNGEN. Let U be a uniformly distri- 
buted, random number and let X be exponentially distributed 
with parameter a. An exponential random number is 
generated by setting u = F(x) and solving for x, 

u = F(x) =1 - e -aX 

. - ax 

1 - u = e 

ln(l-u) = -ax 

x = -ln(l-u)/a, or 

x = -ln(u)/a. 
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A. INTRODUCTION 



MACHREPR is a stochastic, discrete-event simulation of 
a machine-repairman model, which is a special case of 
birth-death processes. The model consists of three heli- 
copters , of which two are in service and one is in cold 
standby, with an option of one or two repairmen. The pro- 
gram output is a graphics display composed of a system 
state-versus -time graph, a table of statistics, and animat- 
ed figures that illustrate the dynamics of the process. 

This manual discusses the model and guides the user in 
the step-by-step operation of the program. In describing 
commands and keyboard entries, the following notation is 
used. Single quotes (' ') enclose commands and phrases 
that are to be entered verbatim. The act of pressing a 
specific key is indicated by brackets (< >) encasing the 
letter or letters on the key. For example, 'format b : ' 
<ENTER> means type the command exactly as it appears 
between the quotes and follow it by pressing the enter key. 
This user's guide shows the commands' in the lower case; 
however, the IBM PC will accept upper case letters as well. 

B. MODEL DESCRIPTION 

MACHREPR contains two versions of the machine-repairman 
model; the one-repairman case and the two-repairman case. 
In both cases, the nominal mission of the model is to keep 
two helicopters flying at all times. If the third 
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helicopter is serviceable, it is kept in cold standby as an 
immediate replacement for the other two. Upon failure of a 
helicopter, a repairman is assigned if he is available; 

otherwise the helicopter joins a queue to await repair. 

Disabled (down) helicopters remain in the repair queue 
until a repairman becomes free. The repair queue is a 
first in, first out type (FIFO). A newly repaired helicop- 
ter will go immediately into flying operation if less than 
two aircraft are in service. If there are already two 
helicopters flying when a repair is completed, then the 
third goes into the cold-standy queue. 

There are several simplifying assumptions for these 

models. All helicopters are assumed to have exponentially 
distributed service lives with the same failure rate a. If 
T is defined as the length of a helicopter's operating 
life, then the mean time to failure (MTTF) is E(T) = l/a. 

The simulation uses the MTTF as the input parameter that 

establishes the particular exponential distribution for the 
operating life. Helicopters will fly continuously until 
failure. A machine that is in cold standby is assumed not 
to be subject to failure. 

Another assumption is that the service time for the 
repairmen has an exponential distribution with parameter p . 
If S is the time required to complete a repair job, then 
the mean time for repair is E(S) = l/p. Like the MTTF, the 
mean time for repair (MTFR) is an input variable for the 
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simulation. The repairmen are assumed to be equally compe- 
tent and thus have the same MTFR. 

In both cases of the machine-repairman model, the state 
of the system is defined as the number of helicopters down 
at time t. The state variable is N(t) = n, n = 0,1,2, 3. 

C. HARDWARE AND SOFTWARE 

1 . Hardware 

The simulation is programmed for the IBM Personal 
Computer (IBM PC) with a color monitor. Since the program 
uses color graphics, the microcomputer must be equipped 
with IBM's Color/Graphics Adapter or a suitable graphics 
card from another manufacturer. MACHREPR will not run 
without the graphics adapter. The program uses two 
assembly-language subroutines that are loaded in memory 
outside of BASIC's 64 kilobyte (K) workspace; therefore, 
there must be a minimum of 128 K of memory installed in the 
microcomputer. The computer must have at least one 
54-inch, floppy disk-drive. 

2 . Software 

MACHREPR and the two subroutines are supplied on a 
distribution diskette. The user must provide the disk 
operating system PC-DOS 2.0 (or 2.10) and Advanced BASIC 
(BASICA) in order to run the simulation. MACHREPR uses 
BASIC's POKE command and certain program transfers, called 
interrupts, to the Basic Input/Output System (BIOS) in the 
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IBM ROM (read-only memory). Because of these interrupts, 
there can be no assurance that the simulation will operate 
properly on other "IBM compatable" machines. 

3 . Program Files 

The distribution diskette contains the following 
program files: 

* MACHREPR . BAS - the main program. 

* RNGEN.SRT - a random-number generator written in 8088 
assembly-language used as a subroutine to the main 
program. 

* SCRNSHFT. SRT - another assembly-language subroutine 
used to perform a specific task on the graphics 
screen . 

* RNGEN.LST - a source-code listing. 

* SCRNSHFT. LST - a source-code listing. 

* MACHREPR . BAT - a batch file that loads and starts the 
simulation . 

The simulation requires the first three files listed above 
in order to run. The two subroutines are binary-image files 
and can not be executed from PC-DOS. The two source-code 
listings are provided as a convenience to the user. 
MACHREPR . BAT is a simple example of how a batch file can be 
used to easily start the simulation. 



D. GETTING STARTED 

1 . Making a Backup Copy 

The first step is to make a copy of the distribu- 
tion diskette. After the copy is made, the original 
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diskette should be stored in a safe place and used only to 
make a replacement for an unserviceable application 
diskette. Although there are several ways that the user may 
choose to create an application diskette, this manual will 
suggest just two: 

* A self -starting (bootable) diskette containing the 
operating system files, the BASICA command file, and 
the simulation's program files. 

* A simple backup copy which contains only the program 
files and which must be used in conjunction with a 
system diskette. 

Owners of microcomputers with a single disk-drive may find 
the first option more convenient. 

To make a bootable diskette, follow the instruc- 
tions in the PC-DOS manual and format a blank diskette 
using the FORMAT command with the /S parameter. The /S 
parameter causes the computer to transfer the operating 
system files from the DOS diskette to the newly formatted 
diskette. The Advanced BASIC file, called BASICA.COM, is 
also found on the DOS diskette and can be copied onto the 
application diskette with the COPY command. Finally, use 
the COPY command again to transfer the simulation's six 
program files from the distribution diskette onto the 
bootable diskette. 

If the second option is desired, then format a 
blank diskette without the /S parameter. This procedure 
does not transfer the operating system files to the appli- 
cation diskette. The next step is to use the DISKCOPY 
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command to make an identical copy of the distribution 
diskette . 

2 . Starting the Simulation 

If the first option was used, insert the bootable 
diskette into drive A and turn on the microcomputer. After 
a few seconds, the DOS prompt A> will appear on the screen. 
It is here that the batch file MACHREPR.BAT can be used. 
Simply type 'machrepr' <ENTER> to start the simulation. 

If using the non-bootable diskette, then turn on 
the microcomputer with the DOS diskette in drive A. Next 
load Advanced BASIC into memory by typing 'basica' <ENTER>. 
After the computer is in BASIC's command mode, remove the 
DOS diskette from drive A and insert the diskette contain- 
ing the simulation. Type 'load "machrepr . bas ", r ' <ENTER> 
to bring the program into memory and to start the 
simulation . 

If desired, the simulation can be started from a 
second disk-drive. Begin as before by turning on the 
microcomputer with the DOS diskette in drive A and loading 
BASICA. With the application diskette inserted into drive 
B, type 'load "b :machrepr . bas " , r ' <ENTER>. During execu- 
tion, MACHREPR reads the program diskette and expects to 
find it in the default disk-drive. When the program fails 
to find the required files on the DOS diskette in drive A, 
it asks the operator for the correct drive label. Simply 
respond to the program's question by pressing <b> . 
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Press any key to continue. . . 



Figure 5: The Title Screen. 



3 . Title Screen 

The display will remain dark for a few seconds 
while the program initializes variables, arrays, and func- 
tion definitions. When initialization is done, MACHREPR 
will present the title screen as shown in Figure 5. 



E. THE PROGRAM MENU 

The program menu, depicted in Figure 6, appears next 
and provides the means of controlling the simulation. With 
the options available on the menu, the user may choose to 
read the program instructions, change the model parameters 
or use the default model parameters. The program instruc- 
tions that are available on the screen are a summarized 
version of those found in this manual. The model 
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parameters that can be altered by the user are the 
helicopter's mean time to failure, the repairmen's mean 
time for repair, and the number of repairmen to employ. 
Additionally, the initial seed for the random-number 
generator can be set. 



PROGRAM MENU 



< I >ns true t i ons . 

<C>hange Model parameters, 

<D>efaul t Model paraMeters, 

<S>et the random numker generator seed, 

Enter your selection. . . 



Figure 6: The Program Menu. 



The simulation starts and ends with the program menu. 
Whenever the user stops the simulation or whenever the sim- 
ulation clock reaches the preset limit of 9950 time units, 
the program returns to the menu. The simulation then can 
be restarted with perhaps a difference set of model parame- 
ters, or the program can be ended altogether. 
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1 . Setting the Random-Number Generator Seed 



To set the seed, press the <s> key. The program 
will then respond by asking for the value of the seed. The 
seed is entered by typing any integer in the range of 1 to 
2,147,483,646 (= 2 31 - 2) and pressing < ENTER > . The pro- 

gram will accept a non-integer value; however, during 
processing the number will be truncated to an integer, with 
the result possibly being different than the desired value. 
After the seed has been entered, the program returns to the 
program menu. 

It is not necessary to set the random-number gener- 
ator seed to run the simulation. There is a default seed 
embedded in the generator's program code, and it is updated 
automatically with each call for a random number. If the 
simulation is stopped and then restarted without setting 
the seed, the generator will use as the initial seed the 
value remaining from the last call in the last simulation 
run . 

2 . Changing the Model Parameters 

When <c> is pressed, MACHREPR will successively ask 
for the desired values of the mean time to failure, mean 
time for repair, and the number of repairmen. The quanti- 
ties are entered by typing a positive integer and pressing 
< ENTER > in response to each question. For best results, 
the values for the mean time to failure and the mean time 
for repair should be at least ten. The program will accept 
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and function properly with non-integer values of MTTF and 
MTFR ; however, these values will be printed on the graphics 
display as integers. Negative numbers will not be accepted 
for the MTTF and the MTFR. The digits 1 and 2 are the only 
valid responses for the number of repairmen. The program 
will immediately start the simulation after the number of 
repairmen has been entered. 

3 . The Default Parameters 

If this option is selected, then the simulation 
starts with the model parameters set at default values. 
These values are: 

MTTF = 15 time units, 

MTFR = 10 time units. 

Number of repairmen = 1. 

The default parameters are selected by pressing <d> . 

4 . Ending the Program 

To end the program, press <e>. MACHREPR will clear 
the screen, exit the graphics mode, and leave the computer 
in the BASIC command mode. If the user wants to exit BASIC 
and return to the DOS, type 'system' <ENTER>. 

F. THE GRAPHICS DISPLAY 

MACHREPR is designed to give the observer the illusion 
of real-time passage as the probabilistic events occur in 
the machine-repairman system. The graphics display indi- 
cates the model's activities. As shown in Figure 7, the 
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graphics display is divided into three sections: the state- 
versus-time graph, the helicopter display, and the table of 
quantities used with queueing models. 




Figure 7: An Example of the Graphics Display. 

1 . The State-Verses-Time Graph 

The upper half of the display contains the state- 
versus-time graph. The values of the state variable N(t), 
the number of machines down at time t, are on the ordinate. 
There are fifty divisions on the abscissa which represent 
units of time. The current time of the system is printed 
at the right-hand end of the abscissa. 

With each unit advancement of the simulation clock, 
the program plots a horizontal line one time unit in length 
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at the height corresponding to the current state of the 
system. The accumulation of these individual plots forms a 



step 


graph as 


shown 


by 


the 


example in 


Figure 


7 . When 


the 


graph 


reaches 


the 


end 


of 


the display. 


the 


program makes 


room 


for the 


next 


plot 


by 


calling SCRNSHFT 


to shift 


the 


entire graph 


one 


unit 


to 


the left. 


In this way. 


the 


model 


's activities 


for 


the 


past fifty 


time 


units can 


be 



seen in a single glance. 

2 . The Helicopter Display 

The three helicopter images in the lower left 
corner of the graphic display indicate the current state of 
the system. The display in Figure 8 shows how the helicop- 
ters are arranged in the initial state 0. All three 
machines are colored green. Helicopters B and C, seen in 
the upper half of the display, are flying. Helicopter A on 
the "ground 1 ' is in cold standby. 

Since all helicopters are initially either flying 
or in cold standby, the next event will be a failure in 
either machine B or C. A machine that fails is grounded 
and transformed into a red colored helicopter with a broken 
rotor. Because there is only one down helicopter, the sys- 
tem is in state 1 and so there is at least one idle 
repairman. The program draws a graphic figure of a man 
beside the down aircraft to represent a machine that is in 
repair. For example. Figure 9 depicts state 1 with 
helicopter B in repair. 
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Figure 8: Initial Helicopter Display in State 0. 




Figure 9: Helicopter Display in State 1. 

As the simulation proceeds, random numbers are gen- 
erated and transformed into interarrival (failure) times 
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and service (repair) times. The program changes the 
appropriate graphic figures when the simulation clock 
reaches the time for the next event. 

Figure 10 shows an example of a helicopter display 
that indicates state 2. In this case, C is in repair. 
Helicopter A is also disabled, and since it does not have a 
repairman, A is in the repair queue. 




Figure 10: Helicopter Display in State 2. 

In Figure 11, all three machines are down, thus the 
system is in state 3. There is a repair queue of length 
two in this example. 

It should be noted that Figures 10 and 11 are 
illustrations from the one-repairman model since only one 
repairman appears in either display. If the simulation is 



48 



running the two-repairman model, then the program will draw 
the second repairman beside the helicopter that fails after 
the one already in repair. 




3 . The Quantities of Interest 

In the study of queueing systems, there are several 
quantities of interest that are used to evaluate the sys- 
tem's performance. In the lower right corner of the 
graphics display (see Figure 7) there is a table of the 
more commonly used quantities. 

The row labels have the following definitions: 

* AR - The average arrival rate of machines into the 

repair system. 

* Wq - The average waiting time of a machine in the 

repair queue. 
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* Lq - The average number of machines in the repair 

queue . 

* W - The average amount of time that a machine is in 

the repair system. This quantity is the sum of 
Wq and the expected service time. 

* L - The average number of machines in the repair 

system. This quantity includes those machines 
in the repair queue and those machines in 
repair . 

* Avl- The machine availability defined as the 

proportion of the time that at least one 
helicopter is flyable. 

* Pd - The probability that a down machine must wait 

for repair. 

* Ps - The probability that an up machine must wait in 

cold standby. 

The numbers that appear in the right-hand column 
under the heading "Limit" are the limiting values of the 
quantities of interest. In other words, the quantities 

converge to the limiting values as time goes to infinity. 
The limiting values are calculated with the balance equa- 
tions that correspond with the particular machine-repairman 
model . 

The values in the left-hand column are estimates of 
the limiting values. These quantities are computed with 
data collected from the simulation. Upon the occurance of 
an event, the program recomputes the estimates after taking 
into account the newly collected data from the last state. 

Because the plots on the s tate-versus -time graph 
must be given in terms of integer coordinates, all time 
variables must be rounded up to the next higher integer. As 
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a result, the estimators are biased since they are based on 
the assumption of continuous time when in fact they are 
calculated with data based on discrete increments of time. 
The effects of the bias are minimal; the simulation still 
provides a sufficient demonstration of the convergence of 
the estimates to the limiting values. 

G. KEYBOARD COMMANDS 

There are three keyboard commands available to the user 
as the simulation is in progress. When <p> (for pause) is 
pressed, the simulation freezes to give the user the oppor- 
tunity to study the graphic display. To continue the simu- 
lation, press <c> . The third command is <s> for stop the 
simulation. The <s> key will send the program back to the 
program menu. 

H. PROGRAM LIMITATIONS 

MACHREPR always starts a simulation under the same con- 
ditions. The initial state is N(0) =0 with helicopters B 
and C in flight and helicopter A in cold standby. There 
are no provisions in the program for starting the 
simulation in another state. Restarting the simulation 
from the program menu causes the current time and all data 
collection variables to be set to zero. 

As noted before, all event times are rounded up to the 
next higher integer for plotting the s ta te-versus -time 
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graph. With exponentially distributed interarrival and 
service times, there is a non-zero probability that the 
sojourn in any particular state will be less than one time 
unit. The program will still alter the helicopter display 
properly and compute the estimates correctly; however, the 
state-versus-time graph may display a jump of two states 
because of the program's inability to plot a line of less 
than one time unit in length. The number of times this 
discrepancy occurs can be minimized if the MTTF and the 
MTFR are set at values greater than ten. 
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APPENDIX B 



MACHREPR . BAS FLOWCHART 



Matn Program 
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Clock Module 
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Clock Module SuDrouttnes 
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Fa f lure Module 
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Part A 
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Failure Module - Part B 
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Failure Module SuOrouttnes 






62 



Repair Module 
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Repair Module - Part A 







STMTKlOt • ♦** 








j 







STDlnCLO* • *0* 







com 


11010 


CTmny c#d holo 
to cola Itonaoy 


1 


! 



STATE • STATE - 1 
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Repair Module - Part 8 
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Repafr Module - Part C 




CALL RHGEMCU) 



PEPEPG TI*e(All) - 
TIME ♦ rSEXPHXKU.HJ) 



rAtLTT>C{At) - TIME 
* FHCXPONTOJ, ALPHA) 




’ 


l 


watthe?** - ** i 






STATE - STATE - 1 



ACM oott tto* 
In Qutut to 
totil »*U tf«# 
In QLj»ue. 



COS Lhfl 131110 


Cn*ng« r 
to roc 
with r# 


“♦a r*lo 
2 r*t o 
04 f r»*n 




WAITfctf 


=«* • ’8* 









CALL G*CEN(U) 



STATE - STATE -1 



RETURN 
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ReDatr Module Subroutines 



SIM! 

SUrouttn# 
13010 










| PUT(X f 14J),PCOKXO* 






1 put(x t io) t cBwno | 







RETURN 



S1«T N 
Subroutine \ 

1M " J 



OoUy routine* 
c«in fro* 1 to 

2 M 







I PUTOT.IO), RTDHfLO | 






| PVKX,Ml),aCD«UH | 







return 
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APPENDIX C 



MACHREPR . BAS PROGRAM LISTING 



0 'MACHREPR. 8AS, Ver 3.1, R. E. Nelsen, 30 Aug 1984 

1 CLEAR : KEY OFF: CLS 

2 BQSU8 30000: 'Check for 8ASICA and color/graphics adaptor 

3 0EF1NT A-Z 

4 DIM P! (3) , SUMP ! (3) ,GRNHEL0 (268) , REDKEL0 (268) ,REDHEL02 (263) 

5 DIM T1MECHRT ( 1609) , HA 1NSCRN ( 8002) 

5 DIM FAILTIME! (3), REPRTIHE! 13) ,START9UE<4) 

7 RNGEM=0: TRUE=- 1 : FALSE-0 : U ! =0: SCRNSHIFT=384 

8 FHT1$=‘II8. tH':FHT2$=*tttt':DRVE$='A: ■ 

9 F 1RSTSCRN=TRUE:HA I TREPR<=*' : FA 1L1NGHEL0$= : STDBYHELQ*=* * 

10 REPA1 REDHEL0i=** : SG5 = STR1 HG« ( 40 , 223) 

50 DEF FNEXP0NT ! (A1 ! ,A2! ) = - L0GIA1II/A2! 

51 DEF FNROUMDUP (A ! ) = 1NT (A!) + 1 

52 DEF FNH1N) (A!,8!,C!)=(-(A!<8! AND A ! <C! > *A ! ) + (- 4B ! < A ! AND B!<C! ) *8! > + (-CC !<A! AND C !<B! )*C! ) 

53 DEF FNHELOPOS1T (A$ ) =- ( INSTRCABC,A*>=l) + (-( 1NSTR('ABC\A*)=2H54> + H 1NSTR t "ABC* , A$) =3) *107) 
55 ' 

50 60SUB 20000: ’Title screen 

80 ON ERROR SOTO 30020: 'Load random number generator. 

82 DEF SEG = iHlAOO: BLOAD DRVEf+ “RNGEN. SRT ’ , 0 

84 ON ERROR GOTO 0 

85 ' 

90 BLOAD DRVE4+ ' SCRNSHFT . SRT * , 38 4 : ‘Load screen shift subroutine 

91 ‘ 

100 '==”=MA1N PROGRAM MENU================================= 

120 SCREEN 0,1: COLOR 14,3,0:W1DTH 40:CLS 

121 LOCATE 7 , 14: PR1 NT 'PROGRAM MENU': PRINT SG* 

122 LOCATE 9,1:PR1NT '(Instructions. ' 

123 LOCATE 1 1 , 1 : PRINT *<C>hange aodel parameters.' 

124 PRINT :PR1NT *<D>efauIt aodel parameters.' 

125 PR1NT:PR1NT *<S>et the random number generator seed.' 

125 PRINT: PRINT *<E>nd the program.': PRINT SG$ 

130 LOCATE 21,2: PRINT 'Enter your selection...'; 

131 ' 

140 G0SU8 10000: AI= INSTRCICDSEicdse'.AI) 

141 IF AI=0 THEN 120 ELSE ON AI GOTO 20100,200,210,150,900,20100,200,210,150,900 

142 ' 

149 ' =====RANDQM NUMBER GENERATOR SEED ROUT 1NE========— ==== 

150 COLOR 14,0 ,0: CLS 

151 LOCATE 1,3: PRINT '» SET RANDOM NUM8ER GENERATOR SEED »': PRINT SG$ 

152 LOCATE 4 ,2 : PRINT 'Permissible seed values are integers' 

153 PRINT ' in the range: 1 to 2147483545.' 

154 LOCATE 7,2: INPUT 'Enter the seed value. ..';A» 

155 IF AK1 OR AD2147483545I THEN LOCATE 7,2: PRINT STRING* (39, 32) : GOTO 154 
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159 DEF SEG = &H1A00 

160 fll ! = INT C Af / 167772 101 >: A2 != INT ( (At-A 1 '* 167772101 > /65536 ! > 

162 'Poke the seed's upper 2 bytes into RN6EN's seed storage. 

163 POKE LH164,A1!: POKE 4H163,A2! 

164 AMAI-A1 H16777210I-A2! *65536 ! :A1 ! = INKAt/256) :A2!= Al-Al! *256 

165 'Poke the seed's lower 2 bytes into RNGEH's seed storage. 

166 POKE & 11162, All: POKE IH161,A2! 

167 GOTO 120: 'Return to aenu 

168 ' 

199 '--PARAMETERS R0UTINE==========================; ===== 

200 COLOR 15,5,0:CLS:HSGI=’>>Enter a positive integer only.<<* 

202 LOCATE 1,6:PR1NT ’* CHANGE MODEL PARAMETERS *’:PRINT SGI 

203 LOCATE 4,1: INPUT ’Enter Mean Tiae To Failure. ..’,A!:1F A!<=0 THEN LOCATE 23,3:PRINT HSGI;:LOCATE 
4,1: PRINT STRINGS (40,32) ; : GOTO 203 ELSE ALPHA!=1/A' 

204 LOCATE 6,1: INPUT ’Enter Mean Tiae For Repair. . . ’,A! : IF A!<=0 THEN LOCATE 23,3:PR1NT HSGI;:LOCATE 
6 , 1 : PRINT STRINGS (40,32) ; :GOTO 204 ELSE HU!=1/A! 

205 LOCATE 23,1:PRINT STRINGS ( 40 ,32) MSBS= ’ ’ 

206 LOCATE 0,1: INPUT ’Enter nuaber of repairaen (1 or 2)...’,NUMRHAN:IF (NUMRHANOl) AND (NUMRHAN02 
) THEN LOCATE 8,1:PRINT STRING! C40.32) ; :G0T0 206 ELSE 250 

210 ALPHA!=1/15:HU!=1/10:NUHRMAN=1 

250 NXTFAIL !=0: NXTREPR ! =0: F IRSTFAILURE=TRUE: NXTEVNT=0: T INE=0:SUMFAIL=0: LASTEVNT=0: STATE=0: XX=31 : SUNW 
0=0 

252 FOR 1=0 TO 3:SUMP! U ) =0: STARTQUE ( I ) =0: NEXT 
260 ’ 

299 '=====THEORETICAL STATISTICS ROUTINE================== 

300 P! (1)=2*ALPHA!/HU! 

320 IF NUMRMAN=1 THEN P! (2)=P! ( 1 ) *P ! (1) :P! (3)=P! (2)*ALPHA!/MU! ELSE P! (2)=P! (1)*ALPHA!/HU!:P! (3)=P! ( 
2 ) *ALPHA ! / (MU ! *MU! ) 

324 P! (0)=1/(1*P!C1)+P!(2>+P!(3>) 

325 P! (1)=P! (1)*P! (0) 

326 P! (2)=P ! (2) *P ! (0) 

327 P! (3)=P! ( 3 ) *P ! (0) 

330 IF NUMRMAN=1 THEN PDOHN!=P! (2)*P! (3):LC!=PD0NN!+P! (3) ELSE LQ!=P! (3):PD0NN'=P! (3) 

340 L!=P!(1)*2»P!(2)+3*P!<3) 

350 LAMBDA ! = (2«P! (0)*2*P! (1)*P! (2) )*ALPHA! 

360 N8!=LQ! /LAMBDA! 

370 W !=L ! /LAMBDA ! 

380 PSTDBV !=P ! (0) 

390 AVL!=P! (0)+(P! (1)+P! (1)*P! (2))/3 
395 ' 

499 ‘=====PRINT MAIN DISPLAY ROUTINE====================== 

500 SCREEN 1,0: COL OR 9,0 

502 IF FIRSTSCRN THEN GDSUB 15000 : F IRSTSCRN=FALSE : GOTO 510 
505 PUT (0,0) ,HA1NSCRN 

510 LOCATE 14 ,27: PRINT USING FHT2S; 1/ALPHA! ; 

511 LOCATE 14,37:PRINT USING FMT2S;l/HUf; 

512 Y=33: 60SUB 10100: ' Print statistics to screen. 
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520 LOCATE 12,37:PR1NT USING FNT2$;T1ME; 

521 GOTO 1000 

522 ' 

899 * "-"PR06RAH END ROUT 1NE=============================== 

900 SCREEN 0,0:W1DTH 80:C0L0ft 4-.PR1NT *Progra» ended. ..‘-.END 

998 - 

999 '=====CLOCK MODULE=============================="===~= 

1000 FAILTIME! 11)=1E+31: 'Calculate 1st failure ti»e (B l C) 

1010 DEF SEG =tH 1 AOO: CALL RNGEN1U!) 

1012 FAILTIME! (2)= FNEIPONT! (U! , ALPHA! ) +TIME 
1020 CALL RNGENfU! ) 

1022 FAILTIME! (3)= FNEXPONT! (U! , ALPHA !) +T1ME 
1030 REPRTIME! ( 1 ) =lE+3 1 : ' 1E*31 = ■infinity 1 
1032 REPRTIME! (2)=1E*31 
1034 REPRTIME ! (3)=1E+31 
1034 STDBYHELO$=*A* 

1099 '** Determine next event *« 

1100 WHILE NXTEVNT<9950: "Sets upper liiit to prg* run ti«e. 

1110 IF NOT STATE=0 THEN 1120 

1112 NXTFA1L != FNMIN! (FAILTIME! (1) , FAILTIME! (2) .FA1LT1HE! (3M 

1114 NXTREPR ! = 1E+31 

1114 GOSUB 10210: GOTO 1180 

1120 IF NOT STATE-3 THEN 1130 

1122 NXTFA1L ! = 1E+31 

1124 NXTREPR!= FNMIN! (REPRTIME! (1) , REPRTIME! (2) .REPRTIME! (3) ) 

1124 GOSUB 1 0220: GOTO 1180 

1130 NXTFA1 L •= FNMIN! (FAILTIME! (1) , FAILTIME 1 (2) .FAILTIME 1 (3) ) 

1132 NXTREPR! 1 FNMIN! (REPRTIME! (1) ,REPRT1ME! (2) .REPRTIME! (3) ) 

1134 GOSUB 10210 
1134 GOSUB 10220 

1180 IF NXTFAIL! <=NXTREPR! THEN 1184 
1182 NXTEVNT = FNRQUNDUP4NXTREPR! ) :GOTO 1200 

1184 NXTEVNT = FNROUNDUP (NXTFAIL ! ) 

1185 ' 

1199 Run clock and update tiie chart ** 

1200 DEF SEG:POKE IH4A,0: 'Clear key buffer 
1210 WHILE TIME < NXTEVNT 

1220 TIHE=T 1ME + 1 

1230 IF T1HE<=50 THEN GOSUB 1 1050: XX=XX+5: FOR 1=1 TO 900: NEXT : GOTO 1240 
1232 DEF SEG = iHlAOO 
1234 CALL SCRNSHIFT 
1238 XX=274:G0SUB 11050 

1240 LOCATE 12,37:PR1NT USING FMT2$;T1ME; 

1241 A$=INKEY$ : IF A$=** THEN 1240 

1242 IF A$= * 5 B OR A$=*S* THEN 120 

1244 IF A$ r, p* OR Ai=*P* THEN GOSUB 10010 
1250 DEF SEG: POKE 4H4A,0: 'Clear key buffer 

1240 WEND 

1241 ' 

1299 ‘tt Update display and determine next event 

1300 SOUND 120 ,3: SUMP ! (STATE ) =SUMP ! (STATE > *T 1ME-LASTEVNT 
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1302 LASTEVNT = NXTEVNT 

1310 IF NXTFfilL! <=NXTREPR! THEN 60SUB 2000 ELSE GOSUB 3000 

1311 ' 

1399 *+ Update statistical estiaates *+ 

1400 LAMBDA !=SUHFAIL/TINE 

1410 IF NUMRMAM= 1 THEN PDOWN ! = T SUMP ! (21+SUMP! (3) ) /TIHE.-LQ! =PDOHN!+SUHP! (31/T1ME ELSE LB!=SUHP! (3)/TI 
HE:PDOHN'=LQ! 

1420 L ! = < SUMP ! (D+SUMP! (2)+SUHP! (21+SUMP! (31+SUMP! (31+SUMP! (311/TIME 

1430 WQ ! =SUMHQ/SU«FAIL : H ! =«Q f +1 /MU ! : PSTDBY !=SUMP ! (0) /TIME 

1440 AVL!=SUHP ! (0) /T1ME+ (SUMP! (D+SUMP ! (1 ) +SUHP! (2) ) / (TIHE+TIHE+TIME) 

1450 Y=25: GOSUB 10100: 'Print estimates to screen 

1451 ' 

1490 A$= INKEY*: IF A$= BB THEN 1500 

1491 IF A*= B s B OR A4- B S B THEN 120 

1493 IF Af=*p B OR AI= B P" THEN 60SUB 10010 

1495 DEF SEG: POKE 4H6A,0: 'Clear key buffer 

1500 HEND : 'No stop conand so go back to line 1100. 

1510 CLS: LOCATE 12,1: PRINT ‘Prograa tiae Iiait reached. 1 
1512 PRINT ‘Press any key to return to aenu... B 

1514 GOSUB 10000: GOTO 120 

1515 ' 

1999 's====fa1LURE M0DULE=============~========— ==—»===== 

2000 SUMFAIL=SUHFAIL+1: STATE = STATE + 1 

2020 IF FIRSTFAILURE THEN 60SUB 1 2100: F IRSTFAILURE=FALSE: GOTO 2130 
2030 AX= INSTR ( B ABC* , FA1L1NGHEL0* ) 

2032 1= FNHELOPOSIT (FA1LINGHEL05 ) 

2040 ON STATE GOTO 2100,2200,2300 

2041 ' 

2100 DEF SEG =4H1A00:CALL RNGENIU!) 

2110 REPRTIME! (AX) = FNEXPONT! (U ! , MU! ) +TIME 
2112 FAILTIME! (AZ)= 1E+31 

2120 6QSUB 12220: '6round flying helo 4 assign repairaan. 

2130 AX- INSTR (* ABC* ,STDBYHELQ*> 

2132 1= FNHELOPOSIT (STDBYHELOJ) 

2140 CALL RNGENIU!) 

2142 FAILTIME! (AX)- FNEXPONT! (U! , ALPHA! )+TIME 

2150 60SUB 12230: ’ Move stand-by helo to flying status. 

2160 STDBYHELO$=" * : RETURN 

2161 ' 

2200 IF NUMRHAN =2 THEN 2240 

2210 GOSUB 12210: "Ground flying helo: no repairaan. 

2220 NAITREPRJ = MID* ("ABC", AX, 1) 

2222 REPRTIME! (AX) = 1 E+31 : FAILTIME ! (AX) = 1E+31 
2224 STARTQUE(AX) = TIME 
2230 RETURN 

2240 GOSUB 12220: ’Ground flying helo 4 assign a repairaan. 

2250 DEF SEG =4H1A00: CALL RNGENIU!) 

2252 REPRTIME! (AX) = FNEXPONT! (U! , HU! )+TIME 
2254 FAILTIME! (AX) = 1E+31: NA1TREPR*= BB 

2260 RETURN 

2261 ' 
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2300 REPRTIHEMAI) = 1E+31 :FAILT IKE ' (AI) * 1E+31 
2310 IF NUHRHAN-2 THEN WAITREPR$= HIDM'ABC'.Al.l) 

2315 STARTQUE (AX) = TIME 

2320 60SUB 12210: 'Ground flying helo: no repairian. 

2330 RETURN 

2331 ' 

2999 '““REPAIR SUBROUT I NE================================ 

3000 AI= INSTR! "ABC* , REPAIREDHELO!) 

3010 * = FNHELOPOSIT (REPAIREDHELO!) 

3020 ON STATE GOTO 3100,3200,3300 

3021 ' 

3100 STDBYHELO! 1 HID* ( ‘ ABC , A2 , 1 ) 

3110 FAILTIHE! (AI) = 1E+31 
3120 REPRTIHEMAI) * IE+31 

3130 GOSUB 13010: 'Change red helo to cold stand-by helo. 

3140 STATE = STATE-1: RETURN 

3141 ' 

3200 DEF SEG =4H1A00:CALL RNGENIU! ) 

3210 FAILTIHE! (AI) = FNEXPONT ! (U ! , ALPHA ! ) +TIHE 
3212 REPRTIHEMAI) = 1E+31 

3220 GOSUB 13020: "Change red helo to green and fly. 

3230 IF NUMRMAN=2 THEN STATE=STATE- 1 : RETURN 
3240 AI- INSTR ( * ABC * , HAITREPR* ) 

3242 I = FNHELOPOSIT (HAITREPR!) 

3244 SUHHQ = SUHHQMIHE-STARTQUE(AI) 

3250 60SUB 13030: 'Change red helo to red helo m/ repairian. 
3260 CALL RNGENIU!): REPRT IKE MAI ) = FNEXPONT! (U! ,HU!)+TIME 
3262 HAITREPR* 1 " 

3270 STATE=STATE-1: RETURN 

3271 ' 

3300 DEF SEG =4H1A00:CALL RNGENIU!) 

3310 FAILTIHE! (AI) = FNEXPONT! IU! , ALPHA! l+TIHE 

3312 REPRTIHE! (AI)= 1E+31 

3320 GOSUB 13020: 'Change to green helo and fly. 

3330 A1I= INSTR I 'ABC*, HAITREPR!) 

3332 X- FNHELOPOS IT (HAITREPR!): 'Helo at head of queue. 

3334 SUHNQ = SUHHQ+T IME-STARTQUE I AI X) 

3340 GOSUB 13030: 'Change red helo to red helo t </ repairian. 
3342 CALL RNGENIU!) 

3344 REPRTIHE! (A1I)= FNEXPONT! (U! ,HU!)+TIHE 

3346 IF NUHRHAN=2 THEN HAITREPR!-":G0T0 3390 

3350 ON AI GOTO 3362,3372,3382 : ' Advance the queue. 

3360 'REPAIREDHELO! = ’A* 

3362 IF HAITREPR!=*B' THEN HAITREPR!=*C’ ELSE HAITREPR!=‘B‘ 

3364 GOTO 3390 

3370 'REPAIREDHELO! = ’B* 

3372 IF HAITREPR!=*A* THEN HAITREPR!='C‘ ELSE HAITREPR!='A' 

3374 GOTO 3390 

3380 ‘REPAIREDHELO! = *C* 

3382 IF HAITREPR!=*A* THEN HAITREPR!-'B* ELSE HAITREPR!*^' 

3390 STATE=STATE-1 : RETURN 
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3391 ' 

9998 ’**m SUBROUTINES *f***M»*mtti*ff»*mtttt*****t*»» 

9999 ' 

10000 DEF SES: POKE LH6A,0: 'General input routine*— ——= 

10001 AS- INKEYS: IF AS*" THEN 10001 ELSE RETURN 

10002 ' 

10010 DEF SEG: POKE LH6A,0: 'Pause routine****— —==—=== 

10011 AS* INKEYS: IF AS=*c’ OR AS*’C’ THEN RETURN ELSE GOTO 10011 

10012 ' 

10099 '^—STATISTICS PRINT SUBROUTINE********** ==«—*« = 

10100 LOCATE 18 , Y: PR 1HT USING FttTli ; LAMBDA ! ; 

10101 LOCATE 1 9 , Y : PR I NT USING FHTISjWQ!; 

10102 LOCATE 20,Y:PRINT USING FMT1S;L8!; 

10103 LOCATE 21,Y:PRINT USING FHT 1 S ; W ! ; 

10104 LOCATE 22,Y:PR1NT USING FHT1S;L! ; 

10105 LOCATE 23, Y:PRINT USING FHT1S;AVL!; 

10106 LOCATE 24,Y:PRINT USING FMT1 $ ; PDOWN ! ; 

10107 LOCATE 25,Y:PRINT USING FHT 1 S ; PSTDBY ! ; 

10108 RETURN 

10109 ' 

10209 '"—NEXT EVENT SUBROUTINES*********** =======— == *== 

10210 IF NXTFAIL i -FAILT IHE ! ( 1 ) THEN FA1LINGHELQ$=*A* ELSE IF NXTFAIL ? =FAILT IHE ! (2) THEN FAIL INGHELOS 
=’B* ELSE FAILINGHELOS=‘C* 

10211 RETURN 

10220 IF NXTREPR‘=REPRT IME ! (1) THEN REPAIREDHELQS= , A’ ELSE IF NX TREPR ! =R£PRT IHE ! (2) THEN REPAIREDHEL 
QS=*B* ELSE REPAIREDHELOS=*C* 

10221 RETURN 

10222 ' 

11049 '=====T1C HARK SUBROUTINE**— ==*=—========**===== 

11050 AX=STATE+1:QH AI GOTO 11051,11052,11053,11054 

11051 Y=75: LINE (XX, Y) - STEP (5,0) ,2: RETURN 

11052 Y=59:LINE (XX, Y) - STEP (5,01 ,2:RETURN 

11053 Y=43:LINE (XX, Yl - STEP (5,0) ,2:RETURN 

11054 Y=27:L INE (XX, Y) - STEP (5,0) ,2:RETURN 

11055 ' 

12099 =====FIRST FAILURE SUBROUT I NE=== ==================== 

12100 DEF SEG =S.H1 AOO: CALL RNGENIU! ) :A!= FNEXPONT 1 (U! ,HU! I+TIHE 

12102 IF FAIL 1NGHEL0S=*B* THEN X=54:X1=78:REPRTIHE! (2) =A! : FAILTIME ! ( 2) =lE+3 1 : GOTO 12106 
12104 IF FAILI NGHELOS=*C* THEN X=107: X 1=131 : REPRT IHE I (3) =A ! : FA1LTIHE ! (3)=1E+31 
12106 PUT (X, 105) , GRNHELO 

12110 PI = 3.141593 

12111 CIRCLE (XI, 161), 12,2,0, PI 

12112 DRAM *C2 BU12 U4 H-12,+10 H-8,-4 DIO BH+20,-16 H+12,+10 R6 rt+4,+5* 

12113 CIRCLE <Xl,165),14,2,7.5tPI/6,I1.5»Fl/6 

12114 DRAW *C2 BD2 L4 NH-2,-15 «-8,-6* : DRAW ’Dll H-4,+5 D4 NL1 NR1 U3 N+6,-3* 

12115 DRAW ’BR20 fl+6,+3 D3 NL1 NR1 U4 H-4,-5 U 1 1 ’ : DRAW ’H-8,+6 NH+2,-15 L4* 
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12116 PAINT STEP (0,4), 2, 2 

12117 DRAM *C0 L4 D6 R1 U6 R6 D6 R1 U6 L4‘ 

12118 CIRCLE STEP (0,4) ,2,0:PSET STEP (0,0), 0 
12120 GET U, 143)- STEP (48,40) ,REDHEL0: XI=Xi+19 

12130 CIRCLE (XI ,165) ,4,1: PAINT (XI, 165), 3,1 

12131 DRAW 'Cl 8)1-4, -2 NL3 M+7,+1 (1-4, -2* 

12132 LINE STEP (0,6)- STEP (2,1),3,BF 

12133 LINE STEP (1,0)- STEP (-4,6) , I ,BF 

12134 LINE STEP (1,0)- STEP (2,6) , I ,BF 

12135 DRAW 'Cl L3 BU9 L2 (1-4, -1 U! M+4,+1 R2‘ : PAINT STEP (-1,1), 1,1 

12140 GET (X, 143)- STEP (4B, 40) , REDHEL02: RETURN 

12141 1 

1220? “"-FAILURE DRAWING SUBR0UT1NES====================== 

12210 PUT (X,105) , GRNHELO: PUT (X,143) , REDHELD: RETURN: ' No repairman 

12211 ' 

12220 PUT IX, 105) , GRNHELO: PUT (X, 143) .REDHELO 

12221 FOR 1=1 TO 300:NEXT:PUT (X, 143) ,REDHELO 

12222 PUT (X, 143) ,REDHEL02: RETURN: 1 With repainan 

12223 ‘ 

12230 PUT (X, 143), GRNHELO: PUT (X , 105) , BRNHELO: RETURN: ' Cold standby to flying status. 

12231 ' 

13009 '——REPAIR MODULE SUBROUT 1NES===================== 

13010 PUT (X, 143) ,REDHEL02: PUT <X,143) , BRNHELO: RETURN : ‘ To cold standby 

13011 ' 

13020 PUT (X, 143) ,REDHEL02: PUT (X, 143) .GRNHELO 

13021 FOR 1=1 TO 200:NEXT: PUT (X, 143) .GRNHELO 

13022 PUT (X, 105), GRNHELO: RETURN: 'Repaired and goes flying. 

13023 ' 

13030 FOR 1=1 TO 200:NEXT:PUT (X,143) .REDHELO: PUT (X, 143) .REDHEL02: RETURN: 'Assign a repairman. 

13031 ' 

14999 ' =====(1AIN SCREEN DRAWING SUBR0UT1NE================== 

15000 LINE <0,8)-(319,8),2:LlNE (0, 10) - (319 , 10) , 2 

15002 LINE (31,23)-(31,80):F0R V=27 TO 75 STEP 16:L1NE (27, Y)- (31 , Y) :NEXT 
15004 LINE (31 , BO) - (28 1,00) 

15006 FOR X=36 TO 281 STEP 5 

15008 IF (X=81) OR (X=131) OR (X=1B1) OR (X=231) OR (X=281) THEN LINE ( X , 80) - ( X , B5) ELSE LINE (X,80) 
-(X,83) 

15010 NEXT 

15012 LINE (0, 100)- (319, 100) ,2: LINE (0, 102)- 1319, 102) ,2 

15014 LINE ( 161 , 102) -(161 , 199) ,2 

15016 LOCATE 1,8:PR1NT * MACHINE - REPAIRMAN MODEL’ 

15018 LOCATE 4,3:PR1NT *3*:L0CATE 6,3: PRINT *2’ :LOCATE 8,3:PR1NT ' 1 *:LOCATE 10,3:PR1NT 'O' 

15020 LOCATE 3, 1:PR1NT *N (t )': LOCATE 12,32:PR1NT 'Ti*e:'; 

15022 ’HelDCopter drawing routine 
15024 PI = 3.!41593:X=25:Y2=165 
15026 G0SU8 1503B 
15028 X=78:Y2=127 
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15030 GOSUB 15038 
15032 1=131 
15034 GOSUB 15038 
15038 GOTO 15082 
15038 Yl=Y2-4 

15040 CIRCLE (I,Y1) ,12,1,0, PI 
15042 DRAW ’Cl BU12 U4 L22 R44' 

15044 CIRCLE U,Y2) ,14,l,7.5*PI/6,11.5»PI/6 
15048 DRAW 'Cl BD2 L4 NM-2,-15 M-8,-6' 

15048 DRAW ‘Dll h-4,+5 D4 NL1 NR1 U3 M+6,-3' 

15050 DRAW 'BR20 tt+8,+3 D3 NL1 NR1 U4 M-4,-5 Ull“ 

15052 DRAW "M-8,+8 NM+2,-15 L4' 

15054 PAINT STEP (0,4), 1,1 

15058 DRAW 'CO L4 D6 R1 U8 R8 D8 R1 U8 LI' 

15058 CIRCLE STEP (0,4) , 1 ,0:PSET STEP (0,0), 0 
15080 RETURN 

15082 'Statistics Routine 

15084 LOCATE 14,22:PRINT "»TTF LOCATE 14,32:PRINT 'MTFR: '; 

15086 LINE (161, 1131-1319, 113), 2 

15068 LOCATE 16,22:PRINT 'Var'-.LOCATE 16,26:PRINT *Esti «ate’ : LOCATE 16,35:PRINT 'Limit' 
15070 LOCATE 18,22:PR1NT 'AR' 

15072 LOCATE 19,22:PRINT 'Wq' 

15074 LOCATE 20,22:PRINT 'Lq' 

15076 LOCATE 21,22:PRINT 'W 
15078 LOCATE 22,22:PRINT 'L' 

15080 LOCATE 23,22:PRINT 'Avl' 

15082 LOCATE 24,4: PRINT “ A* ; : LOCATE 24 , 10:PRINT 'B'; 

15084 LOCATE 24,17:PRINT 'C'pLOCATE 24,22:PRINT 'Pd'; 

15086 LOCATE 25,22:PRIHT 'Ps'; 

15087 GET (0,0) - (319,199) .HAINSCRN 

15088 GET ( 1 , 143) - ( 49, 1 83) ,GRNHELO 

15089 RETURN 

15090 ' 

19999 '=====TITLE SCREEN=========— ======================= 

20000 SCREEN 0,1: COLOR 14,1,1: WIDTH 40: CLS 
20010 LOCATE 2,12: PRINT 'GRAPHIC SIMULATION' 

20012 LOCATE 4,18: PRINT 'of the* 

20014 LOCATE 6,9 : PRINT 'MACHINE - REPAIRMAN MODEL' 

20016 LOCATE 9,14: PRINT 'by R.E. Nelsen' 

20020 LOCATE 14,5: PRINT 'Submitted in partial fulfillment* 

20021 LOCATE 15,2: PRINT "of the requirements for the degree of* 

20022 LOCATE 17,1: PRINT 'Master of Science in Operations Research' 

20023 LOCATE 19,4: PRINT 'from the Naval Postgraduate School' 

20024 LOCATE 20,11: PRINT "Monterey, California" 

20025 LOCATE 22,2: PRINT 'Advisors: J.D. Esary, A.F. Andrus' 

20030 LOCATE 25,2: PRINT 'Press any key to continue. . . " ; : GOSUB 10000 

20040 CLS: LOCATE 1,14: PRINT 'Acknowledgment' 

20041 LOCATE 5,6 : PRINT 'The author would like to thank' 

20042 LOCATE 7 , 2: PR INT 'Associate Professor Bruno Shubert for' 

20043 LOCATE 9,2:PRINT 'the donation of his assembly-language' 
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20044 LOCATE 11,2: PRINT "randoa-nuaber generator which foraed’ 

20045 LOCATE 13,2:PRINT "the basis lor RN6EN.SRT." 

20046 LOCATE 25,2:PRINT "Press any key to continue. . . :60SUB 10000: RETURN 

20047 * 

20099 , ==«=INSTRUCTI0NS«»=“==“™*=««»*~*“==« 

20100 SCREEN 0,1: COLOR 14,1 ,1:WIDTH 40: CLS 

20110 LOCATE 1,4:PRINT "Instructions for the Prograa Menu"; 

20111 LOCATE 3,2: PRINT "Press <s> to set the randoa-nuaber seed" ; : PRINT "and enter a positive intege 
r when asked. The prograa returns to the tenu." 

20112 PRINT: PRINT 1 Press <c> to change aodel paraaeters. The proqraa will ask Tor the aachine's a 
ean tiae to failure (HTTF), the repair-"; 

20113 PRINT "aen's aean tiae for repair (MTFR) , and the nuaber of repairaen to eaploy. The" 

20114 PRINT "HTTF and the HTFR should be positive integers that are greater than 10. The nuaber 

of repairaen can be 1 or 2. The si aul at i on starts after the nuaber of repairaen has been entered. 

■ 

20115 LOCATE 18,2: PRINT "Press <d> to use default paraaeters. "; :PRINT "The siaulation will iaaediat 
ely start with the values: HTTF = 15, HTFR = 10, nuaber of repairaen s 1." 

20116 LOCATE 23,2: PRINT "To end the prograa, press <e>." 

20117 LOCATE 25,2: PRINT "Press any key to continue. .." ; :S0SUB 10000 

20120 CLS: LOCATE 1,13: PRINT "Keyboard Coaaands"; 

20121 LOCATE 4,2: PRINT ’During the siaulation, you aay use the" PRINT "following keyboard coaaands 
, ■ 

20122 LOCATE 7,5: PRINT "<p> - pause the siaulation." 

20123 LOCATE 9,5: PRINT "<c> - continue the siaulation." 

20124 LOCATE 11,5: PRINT "<s> - stop the siaulation and return’ 

20125 LOCATE 13,11: PRINT "to the prograa aenu." 

20126 LOCATE 25,2:PRINT "Press any key to continue. :G0SUB 10000 

20130 CLS : LOCATE 1,13: PRINT "Graphic Display"; 

20131 LOCATE 3,2:PRINT "The variable naaes in the display are" ;: PRINT "defined as:"; 

20132 LOCATE 6,1: PRINT "TIHE: current tiae on siaulation dock."; 

20133 LOCATE 8,1: PRINT ’N(t): no. of aachines down at tiae t . " ; 

20134 LOCATE 10,3: PRINT "AR: average arrival rate of aachines"; 

20135 LOCATE 11,7: PRINT ‘into the repair queue."; 

20136 LOCATE 13,3: PRINT "Hq: average waiting tiae of aachines"; 

20137 LOCATE 14,7: PRINT "in the repair queue.’; 

20138 LOCATE 16,3: PRINT "Lq: average length of repair queue."; 

20139 LOCATE 18,4: PRINT "W; average down tiae of a aachine."; 

20140 LOCATE 20,4: PRINT "L: average no. of down aachines."; 

20141 LOCATE 22,2: PRINT "Avl: aachine availability."; 

20142 LOCATE 25,2: PRINT ‘Press any key to continue...";: 60SUB 10000 

20143 CLS: LOCATE 1,3: PRINT "Pd; probability that a down aachine"; 

20144 LOCATE 2,7: PRINT "aust wait for repair."; 

20145 LOCATE 4,3: PRINT "Ps: probilility that an up aachine"; 

20146 LOCATE 5,7: PRINT "aust wait in cold standby."; 
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20147 LOCATE 8,1: PRINT " The nuabers listed under the coluan heading ‘Lint are the halting va 

lues as tiie goes to infinity of the above defined quantities. These values are coaputed with t 
he balance equations for the aodel." 

20148 LOCATE 15,1: PRINT " The values listed under the coluan heading ‘Estimate' are the estiaate 
s of the baiting values. These quantities are calculated with data collected froa" 

20149 PRINT "the siaulation. The estimates are recoa-puted upon every occurance of an event." 

20150 LOCATE 25,2:PR1NT "Press any key to return to aenu. . . "; : G0SUB 10000:G0T0 120 

20151 ' 

29999 *===s=BASICA AND COLOR/GRAPHICS ADAPTER CHECK========= 

30000 DEF SEG =0: IF (PEEKUH410) AND &H30) <> &H30 THEN DEF SEG:G0T0 30007 

30003 LOCATE 3,1:PR1NT "Sorry..." 

30004 PRINT "You do not have the color/graphics aonitor adapter!" 

30005 PRINT "This siaulation uses graphics and requires that adapter." 

30006 DEF SEG : END 

30007 ON ERROR GOTO 30008: PLAY "pl6":G0T0 30011 
3000B WIDTH 80: CLS: LOCATE 3,1 

30009 PRINT "This siaulation uses advanced BASIC." 

30010 PRINT "Reload this prograa after using the coaaand ’BAS1CA'. " : END 

30011 ON ERROR GOTO 0: RETURN 

30012 ' 

30019 '=====F1LE-L0AD1NG ERROR TRAPPING ROUT 1NE============= 

30020 IF ERR <> 53 THEN ON ERROR GOTO 0 

30022 CLS:LGCATE 3, 1 :PR1HT "On which drive can this prograa be found?" : 60SUB 10000 
30024 DRVE$=A$*" RESUME 82 
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0000 

0000 

0000 

0000 

0001 

0003 

0004 

0006 

0008 

OOOB 

OOOF 

0011 

0013 

0016 

0018 

001A 

001D 

001F 

0021 

0023 

0025 

0027 

0029 



APPENDIX D 



RNGEN. SRT SOURCE-CODE LISTING 



inunuiumiimmunnnnmiitnmnu!} 

RNGEN.SRT - BASIC USR Routine 
by B.0. Shubert 
Modified by R.E. Nelsen 
2 May 1984 

This subroutine generates and returns a unifora 
random-number in BASIC’s single-precision format. 
The seed X is not iaported nor returned to 8ASIC, 
but it is automatically updated. The algorithm 
used is: 

X = AX nod M 



Kith A = 7 A 5 and M = 2 A 31 - 1. RNGEN is called 
froa BASIC using the stateients: 

RNGEN = 0 
DEF SES = JiHlAOO 
CALL RNGEN (U!) 

U! is the returned random nuiber and must be 
defined prior to the first call to RNGEN. 

itsittstsistijissttiisitittssisiissjussiuutsus 



CSEG SEGMENT 

ASSUME CS: CSEG, DS: CSEG, ES: NOTHING 



ORG 

RNGEN PROC 



55 




PUSH 


8B 


EC 


MOV 


IE 




PUSH 


8C 


C8 


MOV 


8E 


D8 


MOV 


fit 


0141 R 


j 

MOV 


88 


14 0143 R 


MOV 


88 


FO 


MOV 


88 


FA 


MOV 


83 


FA 01 


CMP 


77 


4A 


JA 


72 


05 


JB 


3D 


F31D 


CMP 


73 


43 


JNB 


D1 


EO 


HULT32: SAL 


D1 


D2 


RCL 


D1 


EO 


SAL 


D1 


D2 


RCL 


D1 


EO 


SAL 


D1 


D2 


RCL 



0 



FAR 




BP 


;save for BASIC 


BP, SP 


; poi nt tD args on stack 


DS 


; save for BASIC 


AX,CS 


;*ake data accessable 


DS , AX 




AX ,L0SEED 


; get old X 


DX.HISEED 




SI , AX 




DI , DX 




DX, 1 


; check if AIX < M 


MULT 48 


; 48 bit aultiply if not 


MULT32 


;else 32 bit lultiply 


AX , OF31DH 




I1UIT43 




AX, 1 


multiply X by 7 A 5 


DX, i 




AX, 1 




DX, 1 




AX, 1 




DX, 1 
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0026 


2B 


C4 


SUB 


AX, SI 


002D 


IB 


07 


SBB 


OX, 01 


002F 


8B 


FO 


MOV 


SI, AX 


0031 


8B 


FA 


MOV 


01, OX 


0033 


01 


EO 


SAL 


AX, 1 


0035 


01 


02 


RCL 


OX, 1 


0037 


01 


EO 


SAL 


AX, 1 


0039 


01 


02 


RCL 


0X,1 


003B 


01 


EO 


SAL 


AX, 1 


003D 


01 


02 


RCL 


OX, 1 


003F 


2B 


C4 


SUB 


AX, SI 


0041 


IB 


07 


SBB 


OX, 01 


0043 


8B 


FO 


MOV 


SI, AX 


0045 


8B 


FA 


MOV 


01, OX 


0047 


01 


EO 


SAL 


AX, 1 


0049 


D1 


02 


RCL 


OX, 1 


0046 


01 


EO 


SAL 


AX, 1 


004D 


01 


02 


RCL 


OX, 1 


004F 


01 


EO 


SAL 


AX, 1 


0051 


01 


02 


RCL 


DX, 1 


0053 


2B 


C4 


SUB 


AX, SI 


0055 


IB 


07 


SBB 


DX , 01 


0057 


8B 


FO 


MOV 


SI, AX 


0059 


8B 


FA 


MOV 


01, DX 


005B 


01 


EO 


SAL 


AX, 1 


005D 


D1 


02 


RCL 


DX, 1 


005F 


01 


EO 


SAL 


AX, 1 


0041 


D1 


02 


RCL 


DX, 1 


0043 


D1 


EO 


SAL 


AX, 1 


0045 


01 


02 


RCL 


DX, 1 


0047 


26 


C6 


SUB 


AX, SI 


0049 


IB 


07 


SBB 


DX, 01 


004B 


BB 


FO 


MOV 


SI, AX 


004D 


8B 


FA 


MOV 


DI,DX 


004F 


01 


EO 


SAL 


AX, 1 


0071 


01 


02 


RCL 


DX, 1 


0073 


01 


EO 


SAL 


AX, 1 


0075 


01 


02 


RCL 


OX, 1 


0077 


01 


EO 


SAL 


AX, 1 


0079 


01 


02 


RCL 


DX, 1 


007B 


2B C6 


SUB 


AX, SI 


0070 


IB 


07 


SBB 


DX , 01 


007F 


E9 


012F R 


JMP 


UNIF01 


0082 


33 


DB 


i 

MULT48: XOR 


BX, BX 


0084 


8B 


CB 


MOV 


CX,8X 


0084 


01 


EO 


SAL 


AX, 1 


0088 


01 


02 


RCL 


DX, 1 


008A 


01 


03 


RCL 


BX, 1 



; done since Al X^M 
; *ul 1 1 pi y X by 7 A 5 
;by multiplying X by 8 
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; subtract original X 



008C 


DI 


EO 


SAL 


OOBE 


DI 


D2 


RCL 


0090 


DI 


D3 


RCL 


0092 


DI 


EO 


SAL 


0094 


DI 


D2 


RCL 


0096 


DI 


D3 


RCL 


0098 


2B 


C6 


SUB 


009A 


IB 


D7 


SBB 


009C 


IB 


D9 


SBB 


009E 


8B 


FO 


nov 


OOAO 


8B 


FA 


MOV 


00A2 


8B 


CB 


MOV 


00A4 


DI 


EO 


SAL 


00A6 


DI 


D2 


RCL 


00A8 


DI 


D3 


RCL 


00 AA 


DI 


EO 


SAL 


OOAC 


DI 


D2 


RCL 


OOAE 


DI 


D3 


RCL 


OOBO 


DI 


EO 


SAL 


0082 


DI 


D2 


RCL 


0084 


DI 


D3 


RCL 


0086 


2B 


C6 


SUB 


0088 


IB 


D7 


SBB 


OOBA 


IB 


D9 


SBB 


OOBC 


9B 


FO 


nov 


OOBE 


8B 


FA 


nov 


ooco 


8B 


CB 


nov 


OOC2 


DI 


EO 


SAL 


00C4 


DI 


D2 


RCL 


00C6 


DI 


D3 


RCL 


00C8 


DI 


EO 


SAL 


OOCA 


DI 


D2 


RCL 


OOCC 


DI 


D3 


RCL 


OOCE 


DI 


EO 


SAL 


OODO 


DI 


D2 


RCL 


00D2 


DI 


D3 


RCL 


00D4 


2B 


C6 


SUB 


00D6 


IB 


D7 


SBB 


00D8 


IB 


D9 


SBB 


OODA 


8B 


FO 


nov 


OODC 


8B 


FA 


nov 


OODE 


BB 


CB 


nov 


OOEO 


DI 


EO 


SAL 


00E2 


DI 


D2 


RCL 


00E4 


DI 


D3 


RCL 


00E6 


DI 


EO 


SAL 


00E8 


DI 


D2 


RCL 


OOEA 


DI 


D3 


KCL 


OOEC 


DI 


EO 


SAL 



AX, 1 
DX, 1 
9X, 1 
AX, 1 
DX, 1 
BX, 1 
AX, SI 
DX,DI 
BX,CX 

SI, AX ; savi ng new X 

DI, DX 

CX,BX 

AX, 1 ; and repeat this 5 ti«es 

DX, 1 

BX, 1 

AX, 1 

DX, 1 

BX, 1 

AX, 1 

DX, 1 

BX, 1 

AX, SI 

DX , DI 

BX,CX 

SI, AX 

D I , D X 

CX,BX 

AX, 1 

DX, 1 

BX, 1 

AX, 1 

DX, 1 

BX, 1 

AX, 1 

DX, 1 

BX, 1 

AX, SI 

DX,DI 

BX , CX 

SI, AX 

DI, DX 

CX % BX 

AX, 1 

DX, 1 

BX, 1 

AX, 1 

DX, 1 

BX, 1 

AX, 1 
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00EE 


01 02 


RCL 


DX, 1 




00F0 


01 03 


RCL 


BX, 1 




OOF2 


2B C6 


SUB 


AX.SI 




OOF 4 


IB D7 


SBB 


07,01 




00F4 


IB 09 


SBB 


B7.C7 




OOFB 


BB FO 


MOV 


SI, A7 




OOFfi 


SB FA 


MOV 


01,07 




OOFC 


BB CB 


MOV 


C7,B7 




OOFE 


01 EO 


SAL 


AX, 1 




0100 


01 02 


RCL 


07,1 




0102 


01 03 


.RCL 


BX, 1 




0104 


D1 EO 


SAL 


AX, 1 




0104 


D1 02 


RCL 


07,1 




0108 


01 03 


RCL 


B7, 1 




OlOfi 


D1 EO 


SAL 


A7, 1 




010C 


01 D2 


RCL 


DX, 1 




010E 


01 03 


RCL 


BX, 1 




0110 


2B C6 


SUB 


A7.SI 




0112 


IB 07 


SBB 


07,01 




0114 


IB 09 


SBB 


BX, C7 


; BX , DX , AX , now contain AIX 


0114 


D1 E2 


SAL 


07,1 




0118 


01 03 


RCL 


BX, l 


; get K=AIX/2 A 31 to BX 


OllA 


01 EA 


SHR 


D7, 1 


;and X’=A$X(mod 2 A 31) to DX, AX 


one 


03 C3 


ADD 


AX , BX 


calculate X” = X’+K 


01 IE 


S3 02 00 


ADC 


07,0 




0121 


81 FA 7FFF 


CMP 


D7, 7FFFH 


; see if X” < « 


0125 


77 03 


JA 


SUBTRM 




0127 


EB 06 90 


JMP 


UNIF01 


;done if so since X” can't = 


012A 


40 


SUBTRH: INC 


A7 


; el se subtract 11 by adding 1 


012B 


81 DA BOOO 


SBB 


D7.BOOOH 


;and subtracting 2 A 31 


012F 


A3 0161 R 


? 

UNIF01: MOV 


LQSEED, A7 


; stor e new X 


0132 


89 16 0163 R 


MOV 


H I SEED, DX 




0136 


01 EO 


SAL 


AX, 1 


;Divide X by 2 A 31 to aake 


013B 


D1 02 


RCL 


07,1 


;a unifora (0,1). 


013A 


32 C9 


7 OR 


CL, CL 


; Set a counter to zero. 


013C 


01 EO 


FNDEXP: SAL 


AX, 1 


;Shift left and count the 


013E 


01 D2 


RCL 


07,1 


;nuiiber of leading zeros. 


0140 


72 04 


JC 


EXPONT 


;Leading 1 found. 


0142 


FE Cl 


INC 


CL 


; Not found, add 1 to counter 


0144 


EB F6 


JMP 


FNDEXP 


; and look at next bit. 


0144 


01 EA 


EXPONT: SHR 


07,1 


;Shift back to correct value 


0148 


01 OB 


RCR 


AX, 1 


; and suppress leading 1. 


014A 


B3 BO 


MOV 


BL, BOH 


; Compute normalized exponent. 


014C 


2A 09 


SUB 


BL, CL 





BL DX, and AX now contain a single precision, uniform 
(0,1) randof number in BASIC’s floating point foriat. 
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01 4E 


IF 




POP 


DS 


;restore seqaent register 


014F 


8B 7E 06 




NOV 


DI , CBPD+6 


; get addr of variable U 


0152 


8B 25 




NOV 


[Dll, AH 


; pass LSB of iantissa 


0154 


8B 55 01 




NOV 


[ D I ] + 1 , DL 


; pass KSB of aantissa 


0157 


88 75 02 




NOV 


EDI 3+2, DH 


; pass HSB of santissa 


015A 


88 5D 03 




NOV 


EDI 3+3, BL 


;pass exponent 


015D 


58 




POP 


BP 


;restore BP for BASIC 


015E 


CA 0002 




RET 


0 


; FAR return to BASIC 


0161 


4130 


LOSEED 


DU 


4130H 


; storage for l 


0163 


00AB 


HISEED 


DU 


OABH 




0165 




RNGEN 


ENDP 






0165 






CSEB 


ENDS 










END 
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APPENDIX E 



SCRNSHFT . SRT SOURCE-CODE LISTING 



;******************************** t** ********** 

; SCRNSHFT. SRT, Ver 2 - BASIC USR Routine 
;by R.E. Nelsen — 27 July 1984 

; Thi 5 subroutine shifts 5 columns to the left 
; the 4 rows of pixels corresponding to the 
;aodel's state vs tiae graph. It is called 
;froa BASIC using the coaaands: 









; SCRNSHFT = 


0 










; DEF SEB = 


&H1A00 










; CALL SCRNSHFT 










> 

; Thi s subroutines calls the ROM BIOS Type 10H 








; interrupt (Video I/O) using routines 12 i 13. 








;********************************************* 


0000 






CSEB 


SE6MENT 










ASSUME CS:CSEG 


, DS: NOTHING 








I 

; Establ i sh a file header 


for use by the 








; BASIC BLOAD coaaand. 




0000 






HEADER: 






0000 


FD 




DB 


OFDH 


; Code for BLOAD file 


0001 


0000 




DU 


0 


;Seg addr location 


0003 


0000 




DU 


0 


;Qffset location 


0005 


0086 




DU 


RTN_LEN 


Routine length 


0007 






5 

SCRN_SHFT PROC 


FAR 




0007 


55 




PUSH 


BP 


;Save for BASIC 


0008 


FB 




STI 




;Enable interrupts 


0009 


B9 0024 


MOV 


Cl, 36 


; Scr n coluan 36 


oooc 


BA 001B 


A6AIN: MOV 


DX , 27 


;Scrn row 27 


000F 


B4 0D 




MOV 


AH, 13 


;Read pix at row DX 


0011 


CD 10 




1NT 


10H 


;l col CX. Pix in AL 


0013 


83 E9 


05 


SUB 


Cl, 5 


; Shift left 5 col's 


0016 


B4 OC 




MOV 


AH, 12 


; Set = 12 to write 


0018 


CD 10 




INT 


10H 


;pixel to new pos'n 


001A 


83 Cl 


05 


ADD 


C X , 5 


;Shif t right 5 col's 


001D 


83 C2 


10 


ADD 


DX, 16 


; Rom 43 


0020 


B4 OD 




MOV 


AH, 13 


; Set = 13 to read 


0022 


CD 10 




INT 


10H 


;Read pixel 


0024 


83 E9 


05 


SUB 


CX,5 


; Shift left 5 


0027 


B4 OC 




MOV 


AH, 12 


; Set = 12 


0029 


CD 10 




INT 


10H 


;Write pixel 


002B 


83 Cl 


05 


ADD 


Cl, 5 


; Shi f t right 5 col's 



83 



;Now repeat procedure for regaining two rows. 



j 



002E 


83 C2 10 


ADD 


DX, 16 


;Row 59 


0031 


B4 0D 


BOV 


AH, 13 




0033 


CD 10 


INT 


10H 




0035 


83 E9 05 


SUB 


CX, 5 




0038 


B4 0C 


HOV 


AH, 12 




003A 


CD 10 


INT 


10H 




003C 


83 Cl 05 


ADD 


Cl, 5 




003F 


83 C2 10 


ADD 


DL 16 


;Row 75 


0042 


84 0D 


HOV 


AH, 13 




0044 


CD 10 


INT 


10H 




0046 


83 E9 05 


SU8 


a, 5 




0049 


B4 OC 


HOV 


AH, 12 




004B 


CD 10 


INT 


10H 




004D 


83 Cl 06 


ADD 


Cl, 6 


;Shift to new coluin 


0050 


81 F9 01 1 A 


CHP 


Cl, 282 


jls last col reached? 


0054 


72 66 


J8 


A6AIN 


;If no, do next col. 






5 

:Finished with 


rows, so 


now aust blank out last 



;5 coluins of each row. 



i 



0056 


41 


INC 


Cl 


; Col 283 


0057 


B4 OD 


HOV 


AH, 13 


;Set to read dot 


0059 


CD 10 


INT 


10H 


;Read blank dot 


0058 


8A D8 


HOV 


BL,AL 


;Save color code 


005D 


B9 0115 


HOV 


Cl, 277 


; Col 277 


0060 


8A 001B 


NEWCOL : HOV 


DX , 27 


;Ron 27 


0063 


B4 OC 


• HOV 


AH, 12 


; Set to write dot 


0065 


BA C3 


HOV 


AL,6L 


; Set backgrnd color 


0067 


CD 10 


INT 


10H 


; Write blank dot 


0069 


83 C2 10 


ADD 


DX, 16 


;Do again for 


006C 


B4 OC 


HOV 


AH, 12 


;row 43, 


006E 


8A C3 


HOV 


AL,8L 




0070 


CD 10 


INT 


10H 




0072 


83 C2 10 


ADD 


Dl, 16 


;and row 59, 


0075 


B4 OC 


HOV 


AH, 12 




0077 


8A C3 


HOV 


AL,BL 




0079 


CD 10 


INT 


10H 




0078 


83 C2 10 


ADD 


Dl , 1 6 


; and row 75. 


007E 


B4 OC 


HOV 


AH, 12 




0080 


8A C3 


HOV 


AL , BL 




0082 


CD 10 


INT 


10H 




0084 


41 


INC 


Cl 


; New col on right 


0085 


81 F9 01 1A 


CHP 


Cl, 282 


;Last coluin? 


0089 


72 D5 


J8 


NEN_CQL 


; No, so do next col 



;Done, so recover BP end return to SASIC. 



84 



J 



008B 5D 




POP 


BP 


008C CE 




RET 




008D 


SCRNJHFT 


ENDP 




= 0086 


RTNLEN 


EBU 


S - SCRNSHFT 


008D 1A 




DB 


01AH ; Need 


008E 


CSEG 


ENDS 








END 


HEADER 
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